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Abstract It is shown that the interior solution of axially symmetric, station- 
ary and rigidly rotating dust configurations is completely determined by the 
mass density along the axis of rotation. The particularly interesting case of a 
mass density, which is cylindrical symmetric in the interior of the dust con- 
figuration, is presented. Among other things, this proves the non-existence of 
homogeneous dust configurations. 
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1 Introduction 

Dust configurations played an important role in the search for global and phys- 
ically meaningful solutions of Einstein's field equations. Already in the 1920s, 
Lanczos obtained solutions describing rigidly rotating, cylindrically symmetric 
and stationary dust configurations, see [71IS]. A larger class of solutions includ- 
ing the one given by Lanczos was obtained by van Stockum in |12) . These 
solutions describe the interior of all rigidly rotating, axially symmetric and 
stationary dust configurations in terms of an arbitrary solution of a certain 
second order partial differential equation. A closed form of these solutions in- 
volving one arbitrary function is given up to integrations in [6j . Unfortunately, 
this arbitrary function lacks a direct physical interpretation. One intention of 
this paper is to describe the degrees of freedom in the solutions by a physically 
interpretable function, more precisely, the mass density given on the axis of 
rotation. 
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A further generalization was made by Winicour in |14) , where differential 
rotation was considered as well. In this case two functions can be chosen. One 
is completely arbitrary and the other must be an element of the kernel of the 
Laplacian in the flat three dimensional Euclidean space M.^. However, in the 
present paper the attention is turned to rigidly rotating dust, i.e., the van 
Stockum class. 

In Newton's theory of gravity it was shown that isoZaiecfH, axially symmet- 
ric, rotating dust configuration unavoidably collapse to a disk lying in a plane 
perpendicular to the axis of rotation and hence they cannot be stationary, 
see [TD] and references therein. However, this scenario can be prevented by 
distant, stabilizing matter distributions TIT. Therefore, no assumptions about 
the exterior of the dust configuration arc made in the approach presented here. 

In general relativity a similar non-existence theorem is still lacking. But 
some partial results are already known. For instance, in it is shown that 
axially symmetric and stationary dust configurations do not yield an asymp- 
totically flat spacetime provided that the mass density is strictly positive in the 
entire spacetime. This is also the reason, why the singularity in the mass den- 
sity of Bonnor's dust cloud [2], which is a special member of the van Stockum 
class, is inevitable. In this paper this result is generalized to dust configura- 
tions with a boundary seperating it from an arbitrary exterior, but the mass 
density is required to vanish on the boundary, see |15) as well. Furthermore, 
in [1] the non-existence of dust configuration in a spatially compact manifold 
was proven. But note that in both theories, i.e., Newton's theory of gravity 
and general relativity, axially symmetric, stationary and rotating disks of dust 
perpendicular to the axis of rotation exist. An important example is given by 
the Neugebauer-Meinel disk of rigidly rotating dust and its Newtonian limit 
(Maclaurin disks) [9\. 

Bonnor's dust cloud serves as an example providing another interesting 
property of dust configurations in general relativity. The mass density of this 
solution admits a non vanishing gradient along the axis of rotation. This is 
not possible in Newton's theory of gravity, because there the angular velocity 
profile determines the mass density completely. In the case of rigid rotation this 
yields constant mass density. However, we will show that in general relativity 
there is no restriction on the mass density along the axis of rotation. In order 
to do so the Einstein equations for axially symmetric, stationary and rigidly 
rotating dust are solved for an arbitrary real analytic mass density given along 
the axis of rotation. This mass density already yields the interior solution up to 
constants. Hence, the radial profile of the mass density is obtained, too. This 
can be interesting for astrophysical observations. Moreover, it is even shown 
that dust configurations with the Newtonian mass density, i.e., homogeneous, 
do not exist in general relativity. The sole mass density constant along the 
axis of rotation turns out to be the one given by Lanczos in the cylindrical 
symmetric case fT. 



The support of the mass density is assumed to be compact. 
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The paper is organized as follows. Section [5] is devoted to the Newtonian 
case. In particular, attention is paid to the interior solution and its implica- 
tions. It is shown that given the angular velocity curve the mass density is 
uniquely determined and the gravitational potential up to an additive con- 
stant. After the formulation of the problem of rigidly rotating, axially sym- 
metric and stationary dust configurations in general relativity in Sect. [3] some 
non-existence results for dust configurations with a boundary are proven in 
this framework. Afterwards, in Sect. IH the solution of the interior field equa- 
tions is obtained in terms of the mass density on the axis of rotation. These 
results are used in Sect. [S]for a discussion of mass densities, which are constant 
along the axis of rotation. 

Throughout the text geometrical units, in which c — G — 1 holds, are 
chosen. 



2 Dust configurations in Newton's tiieory of gravity 

To compare dust configurations in general relativity with the corresponding 
configurations in Newton's theory of gravity some results regarding stationary, 
axially symmetric rotating dust in the latter theory are recapitulated in this 
section. The matter is characterized by a velocity field v, which vanishes on 
the axis of symmetry, and a mass density n, such that 

v(x,t) = w(p,C)pe^, n{x,t) ^ fi{p,C) (1) 

holds, where w = denotes the angular velocity and {p, V'} are the 

cylindrical coordinates with the corresponding unit vectors {ep,e^,ei^} . The 
gravitational potential U has to satisfy the Poisson equation and the Euler 
equation 

y^U = 47r^, = -AiVC/. (2) 

at 

The field equations and their consequences are investigated in an arbitrary 
open subset Q of the support of the mass density. Hence, it is not necessary 
to exclude other matter distributions in the exterior of J7. Furthermore, the 
behavior at infinity is not restricted and cylindrically symmetric dust configu- 
rations are included in our considerations as well as distant objects stabilizing 
the dust configuration like in jlOj . 

In this section we impose the following smoothness conditions for i? and 
the functions U and fi: 

Assumptions 1 

1. The boundary of Q denoted by df2 is continuously differentiable, 

2. the mass density /x 7^ in fl, 

3. U eC^{[2). 
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Note that the first condition is only for concreteness. Others, not necessarily 
equivalent, like a Ljapunov surface |51 can be considered, too. The third con- 
dition leads with the Poisson equation ([2]) to /i € C"(J7) and thus with the 
Euler equation to cj e C^{f2). Weaker differentiability conditions are possible, 
if weak solutions for the gravitational potential are considered as well. 

From the Eulcr equation ([2]) and the particular form of the velocity field 
([T]) it follows that 

U,^^U.c = 0, U,p^pu{pX)\ (3) 

where a comma denotes a partial derivative. Therefore, the gravitational po- 
tential U is cylindrically symmetric in il and consequently the angular velocity. 
Hence, U is given up to an additive constant in J7 by 

U{p) ^Uo + J puj{pfdp. (4) 

The constant Uq can be fixed using the solution in the exterior of Q and a 
smoothness condition for U across dQ. If ^ is inserted in the Poisson equation 
(HI) all possible mass densities can be obtained given an angular velocity profile 

by 

(w(p)V),p + w(p)' =47rAi(p). (5) 

Therefore, for every angular velocity curve exists one and only one mass den- 
sity. In the case of w = in i7, i.e., rigidly rotating dust, the constant mass 
density 

ZTT 

is obtained. It is worth noting that this proves the non-existence of static dust 
configurations under the assumptions 1 in Newton's theory of gravity. With 
(j4|) the solution is completely determined by the constant wq in fi. 

Conversely, the angular velocity is determined for a given mass density as 
a solution of the ordinary differential equation ([5]) , which is given by 




p{p')p'dp' -f a. (7) 



If the rotational axis intersects i7, the constant a must be chosen to preserve 
the differentiability of the angular velocity in Q. If it does not, additional 
information about the solution in the exterior of J7 is necessary in order to 
determine a. 

The non-existence of isolated, axially symmetric and stationary dust con- 
figurations in vacuum can be shown provided that U G C^(E'^) and U vanishes 
in infinity, see pillO). 
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3 Rigidly rotating dust configurations in general relativity 

In the case of general relativity we restrict ourselves to axially symmetric and 
stationary spacetimes. Hence, it is convenient to use the Lewis-Papapetrou 
line element in quasi cylindrical coordinates 



where the functions U, k, a, W depend only one the coordinates p, C- Since the 
field equations are discussed in a spacetime region Q, where the matter can be 
interpreted as dust, the function W can be chosen to be the radial coordinate 
p by means of a conformal mapping. Furthermore, for rigidly rotating dust a 
transformation in a co-moving coordinate system is possible without changing 
the form of the metric. Let us for simplicity of notation assume the metric 
^ is already given in these co-moving, canonical Weyl-coordinates and let us 
denote with SI an open subset of M^, such that the closure of X? is a subset of 
the restricted coordinate map of Q with respect to {p, ip, C}. 

The only non-vanishing component of the stress-energy tensor in Q reads 



in this coordinate system, where p denotes the non-negative mass density. 
In analogy to the last section we assume the following: 

Assumptions 2 

1. The boundary dfi of SI is continuously differentiable, 

2. the mass density p ^ in S2, 

3. U, a, k& C^{S2). 

Note that again the first condition could be substituted by others like the 
Ljapunov conditions. Furthermore, it is not assumed that the dust configura- 
tion is an isolated object, i.e., that the spacetime is asymptotically flat. Other 
matter distributions can be present in the exterior of Q. In particular, if several 
non-connected components of the dust configuration exist, they can be treated 
independently in the approach to be described. 

If we denote the part of the axis of rotation, which intersects S2, by A, then 
it is convenient for the Theorem [5] to formulate a second set of assumptions: 

Assumptions 3 

1. The set A is not empty, 

2. the origin {p, C.) = (0, 0) lies in A, 

3. the spacetime is elementarily flat. 

If the first condition holds, the second can always be realized by a coordinate 
shift in the ^-direction. If the Assumption I2I3I and the elementary fiatness 
condition are satisfied, then 



ds2 ^ e-2^[e2'=(dp2 + dC^) + W^A^^] - {dt + ad^f, 



(8) 




-2U 



(9) 



a,p, k,p U^p € 0{p) and a_^, k^^ G O(p^) 



(10) 



holds. 
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Now we turn our attention to the field equations. The contracted Bianchi 
identity T"'''.^ = and the Assumptions [5] imply that the function U must be a 
finite constant Uq in f2. Therefore, the non-redundant field equations simplify 
in i7, see, e.g., HI], to 

^{Vaf ^SiTfie^^ (11) 
P P 

where A is the Laplace operator in cylindrical coordinates in the three dimen- 
sional Euclidean space for axially symmetric functions. The function k is given 
by the line integration 



k = 



e4C/o 



1 2 

-((a,c)^ - (a,p)^)dp - -a,pa,<;dC 



(12) 



Equations (jll[) and (TT^ are well defined in a neighborhood of the rotation 
axis because of (fTO)) . 

It is a well known fact that the field equations in the vacuum can be 
simplified to the Ernst equations with the transformation 



5 ^_^e4^, fe.^^e^^. (13) 
P ' P 

The integrability condition of this transformation holds because of the field 
equations (1111) in f2, too, and the function b is twice continuously differentiable 
in f2 because of Assumptions [5] The transformed field equations read 

Ab = 0, (14a) 
(Vfe)^ ^ 87r/ie2''^+2^«, (14b) 

p[((M'-(&c)')dp + 2M,cdC] , (14c) 

where the first equation is the integrability condition of the inverse transfor- 
mation. The behavior of the functions a, k close to the axis (jlOl) ensures that 
the transformation and the field equations (ITil) are also valid on the axis. 
Note that since b is harmonic it is real analytic in f2, as well. With (|14cp 
and (jl4bl) k and /j, are real analytic, too. Therefore, singularities in the mass 
density are excluded by the assumptions and the field equations. Conversely, 
only real analytic mass densities can be given in order to obtain a solution b 
in C'^{f}), which is more restrictive than in Newtonian physics, where only a 
continuous mass density is required. 



4 The solution of the field equations 



Before the general solution is obtained, we prove some non-existence state- 
ments in the formalism presented in the last section. 
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Theorem 1 Let us suppose the Assumvtions \WJ\ and lW^ are satisfied. Then 
the field equations p4p do not admit solutions 6, fc, Uq G C^(i7) H C(i7 U dfi) 
in f2, if one of the following properties is satisfied: 

1. The mass density fi G C^{f2 U df2) vanishes on df2, 

2. b is constant on df2, 

3. the normal derivative of b on dfi vanishes. 

Proof The first case can be reduced to the third using (jl4bp . The cases [2] and 
[3] follow directly from the uniqueness of the solution of the Laplace equation 
for b under the given assumptions. In both cases the unique solution is given 
by 6 = const, which in return leads to = in J7 with the field equation 
(jl4bp and the fact that e^^+'^v is finite and positive, which can be seen from 
the differentiability assumptions. However, this is in contradiction with the 
definition of J7 and the Assumption I2I2I □ 

The first part was also shown in pP. The theorem includes also the fact that 
solutions in the van Stockum class, which describe a spacetime filled completely 
with dust and a mass density vanishing at infinity!^ necessarily have to violate 
some of the assumptions of Theorem [TJ see, e.g,. for other proofs [3 or for a 
recent approach Ij. In particular the differentiability conditions of b are not 
satisfied by the solution given in J^] , because there is a singularity at the origin. 

As we will show in the remainder of this section the conditions which are 
implied by the assumption of rigid rotation are not as restrictive as in the 
case of Newton's theory of gravity. More precisely, to every real analytic mass 
density chosen arbitrarily at A two solutions of the inner field equations in J7 
can be assigned at least locally. Let us denote by the open ball with the 
radius e and the origin of M.^ as center. 

Theorem 2 Let us suppose that Assumptions 2 and 3 hold. Furthermore, let 
us assume that the mass density fx is real analytic in Q in a neighborhood of 
the origin with the radius of convergence e of the series expansion in C. Then 
the solution b G C^(i7) of the field equations (|14p is completely determined in 
a non-empty set d B^ n f2 by the mass density and its derivatives at the 
origin, an arbitrary constant b±{0) and a choice of a sign: 

oo 

b± = 6±(0) ± VSne^Uo ^ -(^)('-i)r'P,(cos6'), (15) 
1=1 

where (y^)*-"'' denotes the nth derivative of with respect to C at the point 
{p,C) = (0,0). The Pi denote the Legendre polynomials of the first kind. Fur- 
thermore, polar coordinates p = rsm9 and C = rcos9 are used. 

^ In order to interpret the last theorem physically the assumptions, which were necessary 
to ensure the validity of the transformation I I13II and which were summarized in the last 
section, have to hold. 
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Proof With and ([T5)) it follows that b p vanishes along A. Hence, the 
second field equation (jl4bp simplifies along A to 

= ±v/87r^e2C^«. (16) 

Therefore, the mass density given along A determines the function 6 up to a 
sign and a constant. Since f2 is an open set and (p, C) — (0, 0) is assumed to be 
an inner point a radius ct > exists, such that <Z B^f] f2. In B^, the square 
root of the mass density at the axis admits the convergent series expansion 

v/^W^Eti^Vm)*'''^'. (17) 

1=0 

Because b is an axially symmetric harmonic function in B„^ it can be written 
in the form 

oo 

6(r,^) =^A,r'PKcos0). (18) 

1=0 

The coefficients Ai can be derived using the identity theorem of power series, 
(HH), (HZl) and (0. This yields the coefficients Ai given in □ 

Some remarks are expedient here. The introduction of the set B^ is for purely 
technical reasons, i.e., to avoid different assumptions, e.g about the topology 
of Q. If the convergence radius e is such that B^ D f2 and fi is a region, then 
the result (ITSl) can be extended to the entire set H. The assumptions about 
the analyticity of fi is necessary and sufficient in order to obtain a solution b 
of the field equations in accordance with the differentiability assumptions. 

To obtain also solutions admitting a non-constant mass density along the 
axis seems surprising in the light of the results in Newtonian gravity in Sect. 
2, where a constant mass density is implied by rigid rotation ([5]). One possible 
explanation is that gravitomagnetic effects due to the motion of the dust will 
act like a force in ^-direction. Thus, other solutions of a generalized Euler 
equation ([2]), at least in a slow motion limit, are possible. These results and 
how to assign to such solutions a proper Newtonian limit using Ehlers frame 
theory will be discussed elsewhere. 

The constants in (jl5p cannot be determined any further. If 6 is a solution 
of the field equations ([T4| in H so are b + const, and —6 and the same mass 
density is obtained from them. The constant b±{0) is the usual freedom due 
to the transformation formulas (jl3p . 

Theorem [2] provides us together with the field equation (jl4cp with an al- 
gorithm to determine the general solution of the field equation if an analytic 
mass density is given along the axis. By (jl5l) 6 is obtained from the mass den- 
sity up to a constant and a sign. With (|14l) the function k as well as the mass 
density can be determined independently of the chosen constant and sign. The 
sign in psp and the constant Uq can be fixed, provided that a solution of the 
field equations in the exterior of f2 is known. 



9 



In order to obtain an exterior solution of the Einstein equations several 
approaches are possible. In some cases the mass density along the axis can 
be extended in ^, e.g., if the radius of convergence of the Taylor expansion in 
C is infinite such that a globally valid cosmological solution can be obtained 
and no exterior solution arises. Such spacetimes describe a universe filled with 
axially symmetric, stationary, rigidly rotating dust. 

If the interior solution described in (fT5|) should be joined to an asymptot- 
ically flat vacuum exterior one has to solve a "Dirichlet problem" with a free 
boundary for the Ernst equation. Whether such a solution exists, especially 
in the light of the non-existence of such dust configurations in the Newtonian 
gravity, is still an open and difficult question. Perhaps our form of the interior 
solution will prove useful to answer it. However, if the vacuum exterior is not 
supposed to be asymptotic flat, then global solutions exist, e.g., van Stockums 
cylindrically symmetric dust [12] or |13j joining the dust to a vacuum exterior. 

Another possibility would be to consider an exterior, where matter can be 
present. A first approach to such an "stabilizing" matter configuration could 
be to consider a shell enclosing the dust configuration. However, one has to 
solve a "Dirichlet problem" for the Ernst equation in the region between the 
dust and the shell and a "Dirichlet problem" for the asymptotic flat vacuum 
region outside the shell. Even though this problem is not trivial the limiting 
case of a shell situated on the surface of the dust seems feasible. These "dust 
stars with a crust" will be investigated in future work. 

5 The non-existence of homogeneous dust configurations 

The algorithm described above is now applied to the important example of 
constant mass density along the rotation axis. It was shown in Sect. [5] that 
this was the sole possible case for rotating Newtonian dust configurations ([3]) 
and that in the case of rigid rotation the mass density must be homogeneous 
([5]) in ]7. As is proven in the following corollary this does not hold in gen- 
eral relativity. The solutions of the field equations in f2 for mass densities 
independent of the C coordinate do not yield a homogeneous mass density. 

Corollary 1 If a dust configuration satisfies Assumptions [D and and the 
mass density ji is constant = ^ along the axis of symmetry, then there 
exists a a > such that the mass density is given in B„ G f2 by 

flip, C) = Moexp(27r^oe-'^V')- (19) 

Proof Because the mass density is constant in A all derivatives with respect to 
C vanish at the origin and the convergence radius e of the series representation 
at this point is infinite. Using the Theorem [5] and the transformation between 
cylindrical and polar coordinates the solution of the field equations (jl4ap and 
(jl4b|) can be written as 



&(p,C) = &±(0)±v/8^e2^«MoC 



(20) 
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in a i?o- C J7 with a > 0. The function k is obtained by means of the hne 
integration (I14cp and dTU]) 

k = -7r/ioe-2^V^ (21) 

Inserting and (|2ip in p4ap yields the mass density given in (IT51) in Bo- . □ 

This coroUary does not only prove the non-existence of homogeneous dust 
configurations, it also gives the only possible mass density in the cylindri- 
cally symmetric case as obtained by Lanczos, see [7 . But here Q need not 
be cylindrically symmetric. Only fi must be independent of C on the axis of 
rotation. 
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